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Abstract 



Microscopic energy current can be de- 
rived from microscopic energy field local- 
ized only in particle positions. If the en- 
ergy current is applied to classical molec- 
ular dynamics, it is expected to become 
a new information. However, multi body 
force except two body force causes obscu- 
rity when multi body interaction energy 
is localized. In the present study, our new 
method enables to determine the localiza- 
tion of the multi body interaction energy. 
We obtain the energy current between 
particles by the law of the multi body 
force corresponding to Newton's third law 
on two body force. 

Keywords ; Energy current, Multi-body force, 
Molecular dynamics 

1 Introduction 

If a physical quantity and its flow exist, they 
need to be fields. The energy should be a field 
if the energy current is considerd. The total 
energy is, however, one quantity in the phase 
space. It is namely not the field. To accept the 
existence of the energy current, microscopic en- 
ergy field is defined which all sort of the energy 
is localized only in particle positions. The en- 
ergy current field is derived from the microscopic 
energy field. S. Lepri, R. Livi and A. Politi de- 
cided the energy current with the approximation 
of low- k limit [I]. S. Takesue derived the energy 
current field which satisfies strictly the 'continu- 
ity equation' with the microscopic energy field 
[2]. In his theory, the energy current field in- 
cludes the term of the energy current between 
particles. The his theory differs from the one of 
J. H. Irving and G. Kirkwood |3] in the point 
that a distribution function is not required. 

The microscopic energy field is not zero only 
at particle positions by the use of Dirac's delta 
function. The kinetic energy of each particle is 
naturally localized in the its particle position. 
Moreover, not only the kinetic energy but also 
the interaction energy is localized in particle po- 



sitions. However, the interaction energy cannot 
be possessed by a single particle because the 
other party is necessary for the interaction. His- 
torically, it was thought that both of two par- 
ticles doing two-body interaction have just half 
the interaction energy. 

If the theory of the microscopic energy field 
is expanded into a classical molecular dynam- 
ics (CMD), the energy current between parti- 
cles can become a new information to analyze 
the dynamics of the particle. When particles 
are atoms and have covalent bonds, they are of- 
ten subjected by the multi-body force in CMD. 
For the energy current, the multi-body inter- 
action energy should be localized in the particle 
positions. The problem of the localization of the 
multi-body interaction energy is not as simple 
as the localization of the two-body interaction 
energy. 

For example, let us consider a water molecule 
composed of two hydrogen atoms and an oxygen 
atom. Three-body interaction is often adopted 
to keep the angle between two covalent bonds 
109.4 degree. The three-body interaction po- 
tential function depends on three atom posi- 
tions. Therefore, the three-body interaction en- 
ergy should be localized in the positions of the 
three atoms. One means of localization is to set 
one-third of the three-body interaction energy 
on each particles. However, we consider this 
means is not correct because the oxygen atom 
and the hydrogen atoms differ in their condition, 
both in the kind of atoms and the magnitude of 
the forces. 

In another scene, such a problem of the lo- 
calization the three-body interaction energy is 
discussed historically ^JElEl- They argue the 
three particles has one-third of the three-body 
interaction energy. However they are ground- 
less. In addition, the problem of the localization 
of the interaction energy has not been discussed 
when it is greater than three-body interaction. 

In the present paper, we propose new method 
of the localization of the multi-body interac- 
tion energy into particle positions. The law of 
the multi-body force corresponding to Newton's 
third law is used. 



The theory of the microscopic energy field 
and the energy current on a three dimensional 
space is introduced in section 2. We derive the 
energy current between particles according to 
only two-body force in section 3. In section 4, 
we describe the new method for the multi-body 
force. And, we discuss our new method in sec- 
tion 5. 

2 Microscopic Energy Field 
and Energy Current 

We introduce the derivation of microscopic en- 
ergy field and microscopic energy current on 
Takesue's work[2j. The Hamiltonian of many- 
particle system is given by 



L — ' Zm 



(2.1) 



where p i and m, are the i— th particle momen- 
tum and mass, U({x}) is the total interaction 
potential which depends on plural particle posi- 
tions {a;}. 

First, we consider the system of particles. 
The i-th particle energy a (t) is localized in the 
i-th particle position Xi. Thereby, the micro- 
scopic energy field e(x,t) is defined by 

e(x,t)=^2ei(t)5(x-Xi(t)), (2.2) 



where S(x) is Dirac's delta function in three di- 
mensional space as 



5{x) = 6(x)S(y)S(z), 



and 



x 




(2.3) 



(2.4) 



The i-th particle kinetic energy p i 2 /2rrii is 
possessed by the i-th particle energy ei(t), and 
is not possessed by efc(i), (k ^ i). Moreover, we 
assume that the total interaction energy U({x}) 
is allocated into particles and is localized only in 
their positions. The interaction energy quota of 
the i-th particle is described by «j. Therefore, 



2nij 



Ui(t). 



(2.5) 



A point to notice is that the interaction en- 
ergy quota Ui(t) is merely an artificial quantity. 
The differential of the interaction energy quota 
Ui(t) cannot generate forces acting on particles. 



The differential of the total interaction poten- 
tial U({x}) corresponds to the force. A main 
problem in the present paper is to determine the 
interaction energy quota Ui(t). We describe the 
method to solve this problem in section 2 and 3. 

Because the space integration of the micro- 
scopic energy field e(x,t) is equal to the total 
energy E in the system, we obtain the following 
condition for the i-th particle energy ej(t), 

E = J dx 3 e(x,t) = ^ (2.6) 

i 

The total energy E obeys energy conservation 
law. 

If microscopic energy current field j(x, t) ex- 
ists, it satisfies the following 'continuity equa- 
tion' with the microscopic energy field e(x,t); 



de(x, t) 
dt 



V-j(x,t) = 0. 



(2.7) 



The time derivative of the e(x,t) becomes 
de(x, t) sr-^ dei(t) 



dt 



dt 



-S(x - Xi(t)) 



^et(t)xi(t)-VS(x-Xi(t))J(2. 



where, in three dimensional space, V5(x) is a 
column vector as 



VS(x) 



( £6(x)S(y)S(z) \ 
6{x)p{y)5{z) 



(2.9) 



When the i-th particle moves with velocity 
Xi, the energy current eiXi is generated and is 
localized only in the i-th particle position Xi 
because the energy is localized only in particle 
positions as Ea. l|2.2|l . Moreover, the energy cur- 
rent between particles is generated by the in- 
teraction between particles. From these point of 
view, we propose the microscopic energy current 
field j(x,t) as follow; 

j(x,t) = y^ei(t)xj(t)S(x - Xi(t)) 

i 

+ 3i->k(i) -j!^-<5(a:; segment, i -> k), 
i,k>i '^ fei ' 

(2.10) 

where ji—,k is the magnitude of the energy cur- 
rent from the i-th particle to the fc-th particle, 
£ ki is relative position vector as 



»ki 



•^k X{. 



(2.11) 



The localization function like the line segment 
from the i-th particle to the fc-th particle is de- 
fined by 



5(x; segment, i — v k) 
= (€ki)x\€ki\$((€ki x (*-aJi))» y 
x S ((€ kl x 



where (• ■ -) y and (• ■ -) z are y- and z-components 
of a vector, and 9 is the Heaviside step function. 
The first term of the right-hand in Eq.J^TOJ 
means the energy current due to the movement 
of particles, and the second term is the energy 
current between particles owing to the interac- 
tion of particles. 

The localization function 5{x] segment, i — > 
fc) has following property; 



V • S(x; segment, i — ► k) 
= 8(x - Xi) - S(x - x k ). 



(2.13) 



From Eas. (|2.1Up and 1|2.13|) . the divergence of 
the microscopic energy current field j(x,t) be- 
comes 

V • j(x, t)=J2 <*(*)*<(*) ' V<J(a: - x,(t)) 

i 

+ J2 ji^k(t)S(x-Xi(t)U2.U) 

where jk^i — —ji^k- The function S(x — Xi(t)) 
appears in second term of Eq. (|2.14|) . because 
the energy current between particle exists only 
on the line segment between particles in (|2.10|l . 
By comparison between Eqs. (|2.8(l . Ij2.14|l . we 
obtain the condition for the energy current be- 
tween particles 



dej(t) 
dt 



5> 

k^i 



k{t). 



(2.15) 



Thus, when the time derivative of the i-th par- 
ticle energy ej(t) can be composed by summa- 
tion V^, we can regard the elements of the sum- 

k^i 

mation as the energy current between particles 
ji—>k- Because the i-th particle energy ei(t) is 
defined by Eq. Ij2.5|) . We must determine the in- 
teraction energy quota Ui (t) to obtain the energy 
current between particles ji^ k ■ 



3 Energy Current with 
Two— Body Force 

We explain the interaction energy quota Ui(t) 
according to two-body interaction. As a re- 
sult, we obtain the energy current between par- 
ticles ji— >fc. The two-body interaction poten- 
tial function is described by 4>ik{rik), where 
fife = \xi — Xk\. We set the following localiza- 



tion with the constants and a ■" t; ; the inter- 
(2.12) action energy localized in i-th particle position 



is 4>ik(rik), and the interaction energy lo- 
calized in fc-th particle position is a k lk ^ik(fik)- 
Thereby, the interaction energy quota Ui was de- 
fined by 



= y^' a i lk) &k(r l k), 



A: 7= 



(3.1) 



where ' means summation about particles 

ky^i 

doing interaction with the i-th one. 

From the right-hand of Eq. (|3.1|1 , the inter- 
action energy quota Ui depends only on canoni- 
cal variables. The kinetic energy p i 2 /2rrii is also 
composed only of canonical variables. There- 
fore, the i-th particle energy ej is the function 
only of canonical variables, the Liouville opera- 
tor can be applied to the time derivative of the 
i-th particle energy ej as 



dt 



= -{H, ei } = -J2{ek,ei}. (3.2) 



By comparison between Eqs. (|2.15|) and (|3.2J) . 

the energy current from the i-th particle to the 
fc-th one ji^k is determined as the element of 
the summation in Eq. I|3.2|) . 

k->k = {e k ,ei} 

Pi du k p k duj 
2rrii dxi 2m k dx k 



(ik)Pi ,Uk)P k 



: (nk) 



dx,. 



.(3.3) 



To explain the constants af k ^ and a k lk ^ , let 
us consider the system of the two particles. The 
difference of the 1-st particle energy e\(t) be- 
tween the times t\ and i 2 is the follows from 
Eqs. ESI and fTT) : 



ei(t 2 ) - ei(h) 



P l (t2? P l (tif 



2r(ii 2nii 
+ at 2) [0i 2 (r 12 (t 2 )) - & 2 (ri 2 (ti))] .(3.4) 

On the other hand, the difference of the 1-st par- 
ticle energy e±(t) between the times t\ and t 2 can 



be described from Eqs. (|2.15l) and (|3.3|) by 



ei(t 2 ) - ei(ti) 



ji^ 2 dt 



(12) (12) \ Pi dcj) 12 {r 12 ) 



a, ' ) • — — — -dt 

J m 

( 12 ) (a. a, \ d(t)i 2 (ri 2 ) 



mi dx\ 



a\ , {x 1 ~x 2 ) 
h d(x 1 -x 2/ 

(12) , (12) \ fPjjh) 2 P^hf 



■dt 



zrrii 2rrii 
+ a[ 12) [<t>i 2 {r X2 {t 2 )) - <Mn 2 (*i))] ■ (3.5) 
For Eqs. I|3.4|l and (|3.5H . the constants satisfy 



(ife) (ife) 
a + ai 



1. 



(3.6) 



Because the both of two particles possess just 
half of the two-body interaction energy histori- 
cally, the constants are determined by 



(ife) (ife) 1 
a l = a k = 2' 



(3.7) 



4 Energy Current with 
Multi-Body Force 

We assume that the interaction energy was lo- 
calized in particle positions even if the multi- 
body force appears. We abandon that the inter- 
action energy quota u, is directly proportionate 
to the multi-body interaction energy as is the 
case with the two-body interaction. We pro- 
pose new method to determine the interaction 
energy quota Ui. 

First, we explain the multi-body interac- 
tion potential. The total interaction potential 
U({x}) is the sum of the multi-body interac- 
tion potential. It can include the two-body in- 
teraction potential. However, it does not include 
the on-site potential. The necessary and suffi- 
cient condition to satisfy the total momentum 
conservation is that the multi-body interaction 
potential depends only on relative position vec- 
tors between particles {£} as 



U({x}) = U({£}). 



(4.1) 



The multi-body force acting on the i— th particle 
due to the variation of the relative position vec- 
tor between the i-th particle and the fc-th one 



£ ik is described by F 
property; 



(ife) 



, and has the following 



dU({t}) ggife 
di lk dx, 



hi 



dxi 



(4.2) 



It is namely that the forces acting on the i-th 
and k-th particles and the one with the varia- 
tion of tik are equal in magnitude and opposite 
in direction. This is the law of the multi-body 
force corresponds to Newton's third law (the law 
of action and reaction). The total multi-body 
force acting on the i-th particle F^ is given by 



dugt}) 

dx. 



E F 

k^i 



(ife) 



(4.3) 



Therefore, the Eq. H4.2JI implies the total mo- 
mentum conservation as 



(ife) 



i k^i 



= EE(*f fe) +*f ) ) 

i k>i 

= 0. (4.4) 

We consider a virtual frame which moves at 
a velocity /3 relative to the original frame as 



Xi(t) = Xi(t) - / f3{t')dt' 



(4.5) 



where quantities in the virtual frame are marked 
tilde as We define the interaction energy quota 
Ui as the negative work of the force acting on the 
i-th particle until time t in the virtual frame by 



Ui{t) 



Fi(t') ■ Xi (t')dt' , (4.6) 



where Fi is the multi-body force acting on the 
i-th particle in the virtual frame as 



dU({t}) 
dsbi 



(4.7) 



Relative position vectors and their deriva- 
tives are invariant under the virtual frame as 



tik tiki 

9t ik d tik 



dx. 



dx. 



= hi - 5, 



(4.8) 
(4.9) 



Because the total interaction potential U({t}) 
depends only on the relative position vectors, it 
is invariant under the virtual frame as 



U({t}) = u({t}). 



(4.10) 



Thereby, the multi-body forces are also invari- 
ant under the virtual frame as 



- (ife) 

F 



dugt}) = dugt}) 

dt ik 9t ik 
^(* k ) \ - ™(ife) _ r 



k^i k^i 



(4.11) 
(4.12) 



For Eqs. I|4.5[) and Q4.12|) . we can describe 
the interaction energy quota Ui of Eq. (|4.6|) by 
the quantities in the original frame as 



u t (t) 



F t (t l )-{x l {t')-f3{t l ))dt'. (4.13) 



where the coefficients Cj and c k are constants. 

By comparison between Eqs. 12.1511 and 
(|4.17(l . the energy current from the i-th parti- 
cle to the fc-th one is determined by 



From Eqs. Ij2.5|l . (|4.13|l and canonical equations 
of motion, the time derivative of the i-th particle 
energy ei become 

^1 = p EL - Fi ■ xi +Ft ■ (3 = Fi • 0.(4.14) 

at rrii 

This implies that the energy current between 
particles ji^k is generated in spite of satisfac- 
tion with the energy conservation law Eq. (|2.6|) 
from Eas. (|4.2jl and [j4.3j) as 



dt 



£e* = £E(jf h) +* , ?°)-0 



i k>i 



= 0. 



(4.15) 



Moreover, the plural virtual frame can be 
employed as long as the energy conservation Eq. 
H2.6JI is satisfied. We apply different velocity (3 ik 
in each force due to the variation of the rela- 
tive position vector £ ik . The interaction energy 
quota Ui(t) is conclusively defined by 

Ul (t) = r^Ff V) . ( Xi (t')-p ik (t'))dt', 

J k^i 

(4.16) 

where Eqs. I|4.5I) and (|4.11|) are used. This new 
definition 1)4. 16[) implies the change of the time 
derivative of the i-th particle energy ei(t) from 
Eq. lj4"Hjl t0 

dd 



dt 



(4.17) 



For the energy conservation law (|2.6|) and Eq. 
(|4.2|l give the following condition to determine 
the velocity /3 ifc ; 



IE-EE* 



(ik) 



(3, 



■i ky^i 



i k>i 

= EEFf ) -(A fc -/3 fci ) 

i k>i 

= 0. (4.18) 

This condition (|4.18(l implies that the velocity 
P ik must be independent of the order of the sub- 
scripts i and k. The velocity f3 ik has the dimen- 
sion of the velocity, and depends on only the i-th 
and fc-th particles. Therefore, we define 



Hik — @ki = c i^i + c kXk, 



(4.19) 



3i 



CkXk) 



(4.20) 



Moreover, when the total interaction poten- 
tial U({£}) is composed of the two-body inter- 
action, the energy current from the i-th particle 
to the k-th one ji~>k <|4.20|) should be consistent 
with Eq. (|3.3(l . Therefore, all of the constants 
Ci and Cfc are determined by 



_ (ik) _ 1 

Cl ~ a k ~ 2 ' 



c k =af k) = l. (4.21) 



The energy current from the i-th particle to the 
fc-th one ji—>k is given by 



1 



Ji^k 



Xi + x k )-Ff. (4.22) 



5 Discussion 

In present work, we extended the theory of the 
microscopic energy field and energy current to 
the system of the multi-body force. The i-th 
particle interaction energy quota Ui is defined 
by Eq. I|4.16(l as the negative work of the forces 
acting on the i-th particle until time t in the 
virtual frames moving at the velocity f3 ik with 
each variation of relative position vector £ ik . As 
a result, the energy current between particles 
ji^k is obtained. 

In the one particle system, the interaction 
potential is the negative work on the particle. 
By analogy with this, at first we defined the in- 
teraction energy quota Ui as the negative work 
of the force acting on the i-th particle until time 
t in the original frame by 



Ui(t) 



J F z (t')-x z (t')dt\ 



(5.1) 



where Fi (t) is the force acting on the i-th parti- 
cle Eq. I|4.3jl . From Eqs. I|2.5(l . (|5.1jl and canon- 
ical equations of motion, the time derivative of 
the i-th particle energy become 



de l 
~dt 



Pi 



El 

mi 



0. 



(5.2) 



Therefore, for the Eq. 1)2. 15[) . the energy current 
between particles ji^k becomes zero. The mi- 
croscopic energy current field j(x,t) is namely 
according to only movement of particles. 

Now, in the macroscopic fluid dynamics, the 
energy current includes a stress term and a heat 



current term besides the term of the movement 
of a fluid particle. If the microscopic energy field 
should be connected with the macroscopic field, 
the energy current between particles ji-tk need 
to be not zero. The definition (|5.1|l is regard as 
a failure. 

We notice that the multi-body interaction 
potential and its forces are invariant under the 
transformation of a reference frame. The virtual 
frame moving at velocity f3 ik relative to the orig- 
inal frame is used. As a result, the energy cur- 
rent between particles is generated. Although 
the time derivative of the i—ih particle energy 
e.i(t) (|4.17l) becomes not zero, the energy con- 
servation law is satisfied owing to the law of the 
multi-body force l|4.2J) corresponding to New- 
ton's third law as Eq. (|4.18jl . 

A problem is the determination of the con- 
stants a[ , ai , Ci and c k . When the interac- 
tion is due to the two-body force, it is suitable to 
split the two-body interaction energy into two 
particles just half as Ea. (|3.7l) . The cause of the 
splitting half is not that two particles has same 
properties like a charge, a mass, a sort of atom, 
and so on. The cause must be that the forces 
on two particles are equal in magnitude due to 
Newton's third law. The constants Ci and c k is 
also determined by 1/2 as Eq. (|4.21() due to the 
law of the multi-body force l|4.2|) corresponding 
to Newton's third law. 

Our method for the multi-body force yields 
new information ji^k to the CMD. However, 
couples of the i-th and A:-th particles for ji^k 
depend on the selection of the independent vari- 
ables {£} of the interaction potential U ({£}). In 
the CMD under the multi-body force, we con- 
sider that the independent variables {£} are cho- 
sen as covalent bonds. 

6 Conclusion 

The theory of the microscopic energy field and 
the energy current was used with only two-body 
force. We propose the new method to deal with 
the multi-body force in the present paper. The 
interaction energy quota Ui in Eq. (|2.5|1 is de- 
fined by Eq. (|4. 1 B|l . It means the negative work 
of the forces acting on the z-th particle until 
time t in the virtual frames moving at the veloc- 
ity (3 ik with each variation of relative position 
vector £ ik . Then, the velocity (3 ik is determined 
by Eq. Ij4.19|l . Our method is based on the law 
of the multi-body force Eq. I|4.2|) corresponding 
to Newton's third law. As a result, the theory of 
the microscopic energy field and the energy cur- 
rent is extended to the multi-body force. The 



energy current between particles ji—,k is derived 
as Eq. 
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